Blackbody radiation shift, multipole polarizabilities, oscillator strengths, lifetimes, 
hyperfine constants, and excitation energies in Ca + 



M. S. Safronova 1 and U. I. Safronova 2 ' 3 
1 Department of Physics and Astronomy, 217 Sharp Lab, 
University of Delaware, Newark, Delaware 19716 
2 Physics Department, University of Nevada, Reno, Nevada 89557, 
3 Institute of Spectroscopy, Russian Academy of Science, Troitsk, Moscow, Russia 

(Dated: November 30, 2010) 

A systematic study of Ca + atomic properties is carried out using high-precision relativistic all- 
order method where all single, double, and partial triple excitations of the Dirac-Fock wave functions 
are included to all orders of perturbation theory. Reduced matrix elements, oscillator strengths, 
transition rates, and lifetimes are determined for the levels up to n = 7. Recommended values 
and estimates of their uncertainties are provided for a large number of electric-dipole transitions. 
Electric-dipole scalar polarizabilities for the 5s, 6s, 7s, 8s, Apj, 5pj, 3dj, and 4d, states and tensor 
polarizabilities for the 4p 3/ / 2 , 5p3/2, 3dj, and 4dj states in Ca + are calculated. Methods are developed 
to accurately treat the contributions from highly-excited states, resulting in significant (factor of 3) 
improvement in accuracy of the 3d 5 / 2 static polarizability value, 31.8(3) a 3 ), in comparison with the 
previous calculation [Arora et al., Phys. Rev. A 76, 064501 (2007)]. The blackbody radiation (BBR) 
shift of the 4s — 3d 5 /2 clock transition in Ca + is calculated to be 0.381(4) Hz at room temperature, 
T — 300 K. Electric-quadrupole 4s — nd and electric-octupole 4s — nf matrix elements are calculated 
to obtain the ground state multipole E2 and E3 static polarizabilities. Excitation energies of the ns, 
np, nd, nf, and ng states with n < 7 in are evaluated and compared with experiment. Recommended 
values are provided for the 7pi/2, 7P3/2, 8pi/2> arm 8P3/2 removal energies for which experimental 
measurements are not available. The hyperfine constants A are determined for the low-lying levels 
up to n = 7. The quadratic Stark effect on hyperfine structure levels of 43 Ca + ground state is 
investigated. These calculations provide recommended values critically evaluated for their accuracy 
for a number of Ca + atomic properties for use in planning and analysis of various experiments as 
well as theoretical modeling. 

PACS numbers: 31.15.ac, 06.30.Ft, 31.15.ap, 31.15.ag 



I. INTRODUCTION 

This work presents a systematic study of atomic prop- 
erties of Ca + ion motivated by its importance for the de- 
velopment of optical frequency standards and quantum 
information processing. 

The current definition of a second in the International 
System of Units (SI) is based on the microwave transi- 
tion between the two hyperfine levels of the ground state 
of 133 Cs. The present relative standard uncertainty of 
Cs microwave frequency standard is around 4 x 10~ 16 
[l[ . More precise frequency standards will open ways to 
more sensitive quantum-based standards for applications 
such as measurements of the fundamental constants and 
testing of physics postulates, inertial navigation, mag- 
netometry, gravity gradiometry, and tracking of deep- 
space probes. Optical frequency standards may achieve 
even smaller relative uncertainties owing to superior res- 
onance line quality factors, allowing shorter averaging 
times and higher stability. Significant recent progress in 
optical spectroscopy and measurement techniques has led 
to the achievement of relative standard uncertainties in 

optical frequency standards that are comparable to the 
Cs microwave benchmark. With extremely low system- 
atic perturbations and better stability and accuracy, such 
optical frequency standards can reach a systematic frac- 
tional uncertainty of the order of 10~ 18 [U, HJ. 



Prospects of optical frequency standard based on the 
metastable 4s — 3^/2 transition in Ca + ion have been 
studied in [J-l6|- In 2009, the first absolute transition 
frequency measurement at the 10 -15 level with a sin- 
gle, lascrcoolcd 40 Ca + ion in a linear Paul trap has been 
reported 0]. The development of an ion clock based 
on Ca + has the technological advantage that all neces- 
sary wavelengths for laser cooling and state manipula- 
tion including lasers for photoionization can be generated 
by commercially available and easy-to-handle solid state 
lasers 0, 0] • The operation of atomic clocks is generally 
carried out at room temperature, whereas the definition 
of the second refers to the clock transition in an atom 
at absolute zero. This implies that the clock transition 
frequency should be corrected in practice for the effect 
of finite temperature of which the leading contributor is 
the blackbody radiation (BBR) shift. Recent experimen- 
tal work Q notes that uncertainty due to BBR shift is 
particularly difficult to improve by experimental means. 
This uncertainty results both from the uncertainty in the 
stability and accuracy of trap temperature measurement 
as well as the uncertainty in the evaluation of the BBR 
shift coefficient (i.e. BBR shift at 300K). In the present 
paper, we improve the accuracy of the BBR shift value 
at 300K by a factor of 3. 

Ca + ions have been used for a number of quantum in- 
formation processing experiments (see Refs. and ref- 
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erences therein) . Most of the elementary building blocks 
for quantum information processing such as state initial- 
ization, long quantum information storage times, univer- 
sal set of quantum logic gates, and readout have been 
realized with high fidelity with trapped ion systems such 
as Ca + [lol - HEj ]. Recently, the merits of a high-fidelity en- 
tangling operation on an optical transition (optical qubit) 
were combined with the long coherence times of two clock 
states in the hyperfine ground state (hyperfine qubit) by 
mapping between these two qubits [9(. Precise under- 
standing of the ion qubit and gate operation decoherence 
properties is aided by precise knowledge of atomic prop- 
erties of this system. 

Properties of Ca + are also of interest to astrophysics 
as the absorption spectrum of the Ca + ion is used to 
explore the structure and properties of interstellar dust 
clouds [HI . Below, we briefly review previous studies 
of Ca + atomic properties. 

The lifetime of the metastable 3d 3 /2 and 3d^/2 levels in 
Ca + was a subject of many theoretical and experimental 
studies owing to importance of these states for the optical 
frequency standards and quantum information applica- 
tions. Early theoretical calculations and measurements 
of the 3 da/2 and 3d§/2 lifetimes in Ca + were reviewed 
in Ref. [l8| . Both high precision measurements and cal- 
culations of the 3dj lifetimes were presented by Kreuter 
et al. in Ref. [18j . A measurement technique was based 
on high-efficiency quantum state detection after coher- 
ent excitation to the 3d 5 m state or incoherent shelving in 
the 3c?3/2 state, and subsequent free, unperturbed spon- 
taneous decay. The calculation of the 3dj — 4s electric- 
quadrupole matrix elements was carried out using an ab 
initio relativistic all-order method which sums infinite 
sets of many-body perturbation theory terms. These 
matrix elements were used to evaluate the 3d radiative 
lifetimes and their ratio [Hj]. In Refs. 0, the rel- 
ativistic coupled-cluster theory was used to perform the 
calculations of these lowest excited 3d 3 /2 and 3d§/2 state 
lifetimes. 

The blackbody radiation (BBR) shift of the 4s — 3d$ 
clock transition, accurate to 3%, and the 3d 5 / 2 tensor 
polarizability were presented by Arora et al. in Ref. [2ll | . 
The calculations were based on the relativistic all-order 
single-double method where all single and double exci- 
tations of the Dirac-Fock wave function are included to 
all orders of perturbation theory. The largest contribu- 
tion to the uncertainty of the BBR shift originated from 
the contribution of the highly-excited nf>? states to the 
3d 5 / 2 static polarizability. 

Relativistic coupled-cluster studies of ionization poten- 
tials, lifetimes, and polarizabilities in singly ionized cal- 
cium was recently presented by Sahoo et al. in Ref. [22| . 
Numerical results were given for the Apj lifetime and the 
4s and 3dj polarizabilities. The polarizabilities of the 4s, 
5s, Ap, bp, 3d, Ad of the Ca + ions were recently evaluated 
by Mitroy and Zhang [23| using the non-relativistic con- 
figuration interaction with semi-empirical core potential 
(CICP) method . 



The first measurement of the 4pz/<>. lifetime in Ca + was 
presented by Smith and Gallagher [24[ using Hanle-effect 
method with optical excitation from the 4s ground state. 
The same technique was used by Gallagher [25| to mea- 
sure the branching ratio of the 4p 3 / 2 decay [17.6(2.0)]. 
The beam-foil technique was used by Andersen et al. (26| 
to measure the oscillator strengths for the 4s — Ap, 
Ap — 5s, Ap — Ad, and Ap — bd transitions. The same 
technique was used by Emmoth et al. [27|- Addition- 
ally, the effects of cascades were analyzed and corrected 
for lifetime measurements. The first pulsed laser excita- 
tion measurements of the Api/ 2 and 4p 3 / 2 level lifetime 
(6.96±0.35 ns and 6.71±0.25 ns, respectively) were re- 
ported by Ansbacher et al. [28]. Gosselin et al. pre- 
sented precision Api/ 2 and 4p 3 / 2 lifetime measurements 
in Refs. [II, [3(1 Two sets of results (6.95±0.18 ns and 
6.87±0.17 ns |i|) and (7.07±0.07 ns and 6.87±0.06 ns 
[30l ]) agree within their uncertainties. The mean lifetimes 
of the 4pi/2 and 4n 3 /2 levels in Ca + were measured by 
Jin and Church [3ll [32J to 0.3% precision using a vari- 
ant of the collinear laser-beam-ion-beam spectroscopy 
technique (7.098±0.020 ns and 6.924±0.019 ns, respec- 
tively). Lifetime of the Ap^o level (6.94±0.18 ns) was 
measured by Rosner et al. [33| using the cascade-photon- 
coincidence technique with a sputtered-atom source. We 
note that recent linearized coupled-cluster calculation 
[2l| disagrees with 0.3% Jin and Church measurement 
by 3%. The accuracy of the calculations should be bet- 
ter than 1% owing to excellent agreement of similar cal- 
culations for all alkali- metal atoms from Li to Fr [34] . 
Mg+dH, Sr+HU, and Ba + [37| with all recent experi- 
mental values. 

Precision measurement of the branching ratios of the 
4^3/2 decay of Ca + was performed recently by Gerritsma 
et al. in Ref. [38j . High precision was achieved by a novel 
technique based on monitoring the population transfer 
when repeatedly pumping the ion between different inter- 
nal states. Authors underlined that forty-fold improve- 
ment was achieved for the A(Ap 3 / 2 — ^ S )/J2j A{Ap 3 / 2 — 
ndj) = 14.31(5) branching ratio in comparison with the 
best previous measurement [38[. 

Warner reported (39[ oscillator strengths for the ns — 
n'p, np — n'd, and nd — n'f transitions with ns = 4s — 9s, 
np = Ap — 9p, nd = 3d — Id, and nf = Af — If. The ra- 
dial wave functions were calculated using scaled Thomas- 
Fermi-Dirac wave functions with including spin-orbit in- 
teraction [39(. The relativistic pseudopotential approach 
was applied by Hafner and Schwarzt [4y] to the calcula- 
tion of the ns — n'p electric-dipole transition probabilities 
for n = A — 7 and n' = 4 — 6. Semiemerical wave func- 
tions were used by Theodosiou [4l[ to evaluate oscillator 
strengths and lifetimes of the 5s, 6s, 7s, Ap, bp, 6p, Ad, bd, 
and Af levels. Relativistic many-body theory was applied 
by Guet and Johnson [42| to determine amplitudes for the 
Ap — As and Ap — 3d transitions and the Ap lifetime. Multi- 
configuration Hartree-Fock wave functions were used by 
Vaeck et al. [43|] to evaluate oscillator strengths for the 
Ap — As and Ap — 3d transitions and the Ap lifetime in 
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Ca + . Liaw [44J performed ab initio calculations based on 
the Brueckner approximation for the amplitudes of the 
4p — 4s and 4p — 3d transitions and the 4p lifetime in 
Ca+. The lifetimes of the 5s, 6s, 4p, 5p, 6p, 4d, 5d, and 
4/ levels in Ca + were evaluated by Melendez et al. [45[ 
using the Thomas-Fermi-Dirac central potential method 
in the frozen core approximation and including the polar- 
ization interaction between the valence electron and the 
core using a model potential. Recently, non-relativistic 
CICP method was used by Mitroy et al. |46j to evaluate 
the 4s — np, Ap — ns, 3d — np oscillator strengths with 
n = 4, 5. 

The hyperfine structure of the 4s, 4p, and 3d states 
was evaluated by Martensson-Pcndrill and Salomonson 
using many-body perturbation theory. Numeri- 
cal values were given for the A(4s), A(4p), and A(3d) 
magnetic-dipole hyperfine constants and for the B(Ap^/o) 
and B(3d) electric-quadrupole hyperfine constants [471 ]. 
First measurements published several years later [48| con- 
firmed theoretical predictions [47j . The experimental and 
theoretical results for the A(4s 1/2) and A{Apj) hyperfine 
constants agreed at the 3% level [HI]. Additional cor- 
relation contributions were added to the method used 
in [47j by Martensson-Pendrill et al. [49| to recalcu- 
late all of the above magnetic-dipole hyperfine constants. 
The hyperfine structure splittings of the 4s, Api/2, and 
4p3/2 levels in 43 Ca II were measured by fast ion beam 
collinear laser spectroscopy in (Hoj . Precise determina- 
tion of the 4s ground state hyperfine structure splitting 
of 43 Ca+ (3225.6082864(3)MHz) was reported by Arbes 
et al. [51]. The Doppler-free and potentially very nar- 
row resonances were used in Ref. [52j to determine the 
magnetic dipole hyperfine interaction constant A for the 
Apx/2 and 3d 3 / 2 states of 43 Ca + . Hyperfine structure in 
the three 3d 3 / 2 .5/2 =^ 4pi/2,3/2 transitions were studied 
by fast ion beam collinear laser spectroscopy for all sta- 
ble Ca isotopes in Hyperfine structure parameters 
A(Ap), A(3d), B(4p 3 / 2 ), and B(3d j ) for the odd isotope 
43 Ca + , as evaluated from the splittings observed, agreed 
well with theoretical predictions from relativistic many- 
body perturbation theory [53, Ei| . Recently, relativistic 
many-body calculations were performed to calculate the 
magnetic-dipole hyperfine constants A(As), A{Ap), and 
A(3d) and the electric quadrupole constants -B(4p 3 / 2 ) 
and B{3d) in 43 Ca II 0, The relativistic coupled 
cluster theory was employed by Sahoo et al. 54] to cal- 
culate the hyperfine A constants of the low-lying states. 

Despite many previous studies, no reliable recom- 
mended values exist for a number of properties of low- 
lying states of Ca + . In many cases, only semi-empirical 
calculations are available. In this work, we carry out 
a systematic study of Ca + energies, El, E2, E3 ma- 
trix elements, transition rates, lifetimes, A and B hyper- 
fine constants, E2 and E3 ground state polarizabilitics, 
scalar El polarizabilities of the 5s, 6s, 7s, 8s, Apj, 5pj, 
3dj, 4dj states, and tensor polarizabilities of the 4p 3 / 2 , 
5^3/2; 3d, , and Adj states using a high-precision all-order 
approach [57| in several different approximations. We 



evaluate the uncertainties of our calculations for most of 
the values listed in this work. We also re-evaluated the 
blackbody-radiation shift in a Ca + ion optical frequency 
standard and improved its accuracy by a factor of 3. The 
quadratic Stark effect on the hyperfine structure levels of 
43 Ca + ground state is investigated. The methodologies 
for evaluating the uncertainties of theoretical values cal- 
culated in the framework of the all-order approach are 
developed. The calculation of uncertainties involved es- 
timation of missing high-order effects and ab initio cal- 
culations in different approximations to establish the size 
of the higher-order corrections and approximate missing 
contributions. We evaluated the uncertainties of the rec- 
ommended values for the transition matrix elements, os- 
cillator strengths, transition rates, lifetimes, polarizabil- 
ities, BBR shift, and the Stark shift coefficient. 

The main motivation for this work is to provide rec- 
ommended values critically evaluated for their accuracy 
for a number of atomic properties via a systematic high- 
precision study for use in planning and analysis of various 
experiments as well as theoretical modeling. 



II. ENERGY LEVELS 

Energies of the nlj states in Ca II are evaluated for 
n < 7 and I < 3 using both third-order relativis- 
tic many-body perturbation theory (RMBPT) and the 
single-double (SD) all-order method. The all-order (lin- 
earized coupled-cluster) method and its applications are 
discussed in detail in review [5 71 ] and references therein. 
Therefore, we do no repeat the method description in this 
work, with the exception of the details needed to discuss 
the evaluation of uncertainties in the transition matrix 
elements carried out in the next section. The compari- 
son of the calculated energy levels with the experimental 
values gives an excellent indication of the accuracy of 
the approach and identifies state with particularly large 
correlation corrections. Moreover, a number of highly- 
excited energy levels that we need for the polarizabil- 
ity calculations appear to be not known and are missing 
from the NIST database [56j. Results of our energy cal- 
culations are summarized in Table |TJ Columns 2-7 of 
Table [J give the lowest-order DF energies E^ , second- 
order and third-order Coulomb correlation energies 
and first-order and second-order Breit corrections 

£?W and B^ 2 \ and an estimated Lamb shift contribu- 
tion, E( LS \ The Lamb shift i?( LS ) is calculated as the 
sum of the one-electron self energy and the first-order 
vacuum-polarization energy. The vacuum-polarization 
contribution is calculated from the Uehling potential us- 
ing the results of Fullerton and Rinker (58[ . The self- 
energy contribution is estimated for the s, px/2 an d p 3 /2 
orbitals by interpolating among the values obtained by 
Mohr [H, HH using Coulomb wave functions. For 
this purpose, an effective nuclear charge Z e g is obtained 
by finding the value of Z c ^ required to give a Coulomb 
orbital with the same average (r) as the DF orbital. We 
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TABLE I: Zeroth-order (DF), second-, and third-order Coulomb correlation energies E^ n > , single-double Coulomb energies 



E , ^extrai first-order Breit and second-order Coulomb-Breit corrections B^ n > to the energies of Ca II. The total energies 
(£ ( t 3 't = E<M + £ (2) + E^ + + + EV*\ E?° = E^ + E SD + E^ la + + B ^ + £< LS >) of Ca II are compared 
with experimental energies -Enist HH, SE = E to t - -Enist- Units: cm -1 . a Recommended values. 
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2 


7*1/2 


-16053 


-288.7 


54.1 


1.7 


-1.6 


0.0 


-16287 


-268.8 


32.1 


-16289 


-16304 


16 


14 


7P1/2 


-14051 


-221.0 


35.8 


1.9 


-1.2 


0.0 


-14236 


-218.3 


22.8 


-14246 


























- 14253° 








7P3/2 


-14032 


-219.2 


35.5 


1.4 


-1.3 


0.0 


-14215 


-216.6 


22.6 


-14225 


























-14232° 








7rf 3 /2 


-10576 


-278.3 


46.3 


1.0 


-2.6 


0.0 


-10809 


-250.2 


27.3 


-10799 


-10818 


9 


18 


7^5/2 


-10573 


-277.6 


46.2 


0.6 


-2.6 


0.0 


-10806 


-249.7 


27.2 


-10796 


-10815 


9 


18 


7/5/2 


-8974 


-51.8 


5.5 


0.0 


0.0 


0.0 


-9020 


-54.9 


5.2 


-9023 


-9025 


5 


2 


7/7/2 


-8974 


-51.8 


5.5 


0.0 


0.0 


0.0 


-9020 


-54.8 


5.2 


-9023 


-9025 


5 
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find that the values of E^ LS ^ are very small. For states 
with / > 0, the Lamb-shift is estimated to be smaller than 
0.1 cm -1 using scaled Coulomb values and is negligible 
at the present level of accuracy. We list the all-order 
SD energies in the column labeled E SB and the part of 
the third-order energies missing from i? SD in the column 
labeled ^xtra- The sum of the seven terms E^°\ _E SD , 
E^ tw and gives our final all-order re- 

sult Ef® , listed in the eleventh column of Table Q] Rec- 
ommended energies from the National Institute of Stan- 
dards and Technology (NIST) database [56[ are given in 
the column labeled -Enist ■ Differences between our third- 
order and all-order calculations and experimental data, 
5E<® = E^l - E msT and SE SD = - £ NIST , are 
given in the two final columns of Table HI respectively. 

As expected, the largest correlation contribution to the 
valence energy comes from the second-order term E^' 2 \ 
Therefore, we calculate E^ with higher numerical ac- 
curacy. The second-order energy includes partial waves 



up to l maK = 8 and is extrapolated to account for con- 
tributions from higher partial waves (see, for example, 
Refs. [62|,l63| for details of the extrapolation procedure). 
As an example of the convergence of E^ with the num- 
ber of partial waves I, consider the ground 4s state. Cal- 
culations of £ (2) with ; max = 6 and 8 yield E^ (4s) = 
—4726.2 and —4743.7 cm" 1 , respectively. Extrapolation 
of these calculations yields —4786.3 and —4786.9 cm -1 , 
respectively. Therefore, the numerical uncertainty in the 
second-order value Z?( 2 )(4s) is 0.6 cm -1 . It should be 
noted that the 17.5 cm -1 contribution from partial waves 
with I > 6 for the 4s state is the largest among all states 
considered in Table HJ smaller (about 4 — 6 cm" 1 ) con- 
tributions are obtained for the 3d, 4p, and 4c? states and 
much smaller contributions (0.5 — 1.5 cm" 1 ) are obtained 
for the n = 6 states. 

Owing to complexity of the all-order calculations, we 
restrict I < ^ nax = 6 in the E su calculation. The second- 
order contribution dominates E SD ; therefore, we can use 
the extrapolated value of the E^ described above to 



5 



account for the contributions of the higher partial waves. 
The partial waves I < 6 are also used in the calculation of 
Since the asymptotic I- dependence of the second- 
and third-order energies are similar (both fall off as 
we use the second-order remainder as a guide to estimate 
the remainder in the third-order contribution. The term 
.©extra m Table HI which accounts for the part of the third- 
order MBPT energy missing from the SD expression for 
the energy, is smaller than by an order of magnitude 
for the states considered here. 

The column labeled SE SD in Table U gives differences 
between our ab initio results and the experimental values 
. The SD results agree significantly better with mea- 
sured values than do the third-order MBPT results (the 
ratio of 5E^ /5E SB is about 10 for some of cases), illus- 
trating the importance of fourth and higher-order corre- 
lation corrections. 

We provide recommended values for the lpi/2 arid 
7^3/2 energies in Table U in the separate rows. We es- 
timate these values to be accurate to about 3 cnr 1 . Our 
recommended values for the 8pi/ 2 and 8p 3 / 2 energies are 
— 10221 cm -1 and —10209 cm -1 , respectively. 



III. ELECTRIC-DIPOLE MATRIX ELEMENTS, 
OSCILLATOR STRENGTHS, TRANSITION 
RATES, AND LIFETIMES IN CA II 



these terms). 

Next four columns contain four different all-order cal- 
culations. Ab initio electric-dipole matrix elements eval- 
uated in the all-order SD (single-double) and SDpT ap- 
proximations (single-double all-order method including 
partial triple excitations [34j ) are given in columns la- 
beled Z SD and Z SE) p t of TablcE] The SD and SDpT ma- 
trix elements Z srj include completely, along with im- 
portant fourth- and higher-order corrections. The fourth- 
order corrections omitted from the SD matrix elements 
were discussed by Derevianko and Emmons [65| . Differ- 
ence between the Z su and Z SDpT values is about 0.5 % 
- 2.0 %. 

We have developed some general criteria to establish 
the final values for all transitions and evaluate uncertain- 
ties owing to the need to analyze a very large number of 
transitions. To evaluate the uncertainties of our matrix 
element values and to provide recommended values, we 
carried out semi-empirical evaluation of the missing cor- 
relation corrections using the scaling procedure described 
below. 

The matrix elements of any one-body operator Z = 
y\„. Zij a]a,j are obtained within the framework of the 



SD all-order method as 



(#u,|£|*„) 



y/{* v \*„)(* w \* w y 



(1) 



A. Electric-dipole matrix elements 



where and \^ w ) are given by the expansion 



In Table HH we list our recommended values for 58 El 
ns — n'p and nd — n'p transitions. We note that we have 
calculated over 500 El matrix elements to evaluate polar- 
izabilities and BBR shift presented in this work. We refer 
to these values as "best set" of the matrix elements. We 
list only the matrix elements that give significant contri- 
butions to the atomic properties calculated in the other 
sections. To evaluate the uncertainties of these values, we 
carried out several calculations in different approxima- 
tions. To demonstrate the size of the second, third, and 
higher-order correlation corrections, we list the lowest- 
order Dirac-Fock (DF) Z DF , second-order Z( DF+2 \ and 
third-order ^( DF + 2 + 3 ) values in the first three numeri- 
cal columns of Table [TTJ The absolute values in atomic 
units (aoe) are given in all cases. The many-body per- 
turbation theory (MBPT) calculations are carried out 
following the method described in Ref. 64]. The val- 
ues Z( DF + 2 ) are obtained as the sum of the second-order 
correlation correction Z^ 1 and the DF matrix elements 
Z DF . The second-order Breit corrections B^ are very 
small in comparison with the second-order Coulomb cor- 
rections Z^ (the ratio of B^ to Z^ are about l%-2%). 
The third-order matrix elements z( DF+2+3 ) include the 
DF values, the second-order Z^ results, and the third- 
order Z^ correlation correction. Z^ includes random- 
phase-approximation terms (RPA) iterated to all orders, 
Brueckner orbital (BO) corrections, the structural radi- 
ation, and normalization terms (see [HJ for definition of 



1 ~T~ ^ ^ Prna^jn^a cy ^ ^ Pmnab^tr, 



ala b a a 



mnab 



,(°)\ 



Et t 
Pranva ^rn^rS^ a 



l*! 0) >, (2) 



and ) is the lowest-order atomic state vector. In 
Eq. @, the indices m and n range over all possible vir- 
tual states while indices a and b range over all occupied 
core states. The quantities p ma , p m v are single-excitation 
coefficients for core and valence electrons and p m nab and 
Pmnva are double-excitation coefficients for core and va- 
lence electrons, respectively. In the SD approximation, 
the resulting expression for the numerator of Eq. (fTJ) con- 
sists of the sum of the DF matrix element z wv and 20 
other terms that are linear or quadratic functions of the 
excitation coefficients. The all-order method yielded re- 
sults for the properties of alkali-metal atoms and many 
other monovalent systems [3414371 . |57| in excellent agree- 
ment with experiment. For example, the SD results for 
the primary ns — npj El matrix elements of alkali-metal 
atoms agree with experiment to 0.1%-0.5% [34J. How- 
ever, triple corrections are important for many of the 
nd — n'p matrix elements and have to be included. Our 
ab initio SDpT values include corrections to the equa- 
tions for the valence excitation coefficients p mv and va- 
lence energy. These corrections arise from the addition of 
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TABLE II: Recommended values of the reduced electric-dipole matrix elements in atomic units. The first-order, second-order, 
third-order MBPT, and all-order SD and SDpT values are listed; the label "sc" indicates the scaled values. Final recommended 
values and their uncertainties are given in the Z final column. The last column gives relative uncertainties of the final values in 
%. Absolute values are given. 



Transition 


Z DF 


^(DF+2) 


^(DF+2+3) 




7 (SD) 


^SDpT 


^SDpT 


^■final 


Unc. (%) 


4si /2 


4fl/2 


3.2012 


3.0045 


2.8826 


2.8978 


2.9071 


2.9131 


2.9071 


2.898(13) 


0.45 


4S1/2 


4p 3 /2 


4.5269 


4.2499 


4.0773 


4.0989 


4.1119 


4.1204 


4.1119 


4.099(18) 


0.45 


5Si/ 2 


4pi/ 2 


2.1084 


2.1458 


2.0799 


2.0660 


2.0735 


2.0706 


2.0717 


2.073(11) 


0.51 


5Si/ 2 


4p 3 / 2 


3.0142 


3.0653 


2.9752 


2.9551 


2.9647 


2.9614 


2.9622 


2.965(14) 


0.46 


5S1/2 


5P1/2 


6.4426 


6.3777 


6.1965 


6.2195 


6.2297 


6.2392 


6.2287 


6.23(1) 


0.16 


5S1/2 


5P3/2 


9.1006 


9.0095 


8.7523 


8.7850 


8.7999 


8.8130 


8.7984 


8.80(1) 


0.17 


6S1/2 


4P1/2 


0.5798 


0.6002 


0.5837 


0.5817 


0.5820 


0.5825 


0.5827 


0.582(1) 


0.18 


6S1/2 


4p 3 / 2 


0.8239 


0.8522 


0.8291 


0.8264 


0.8267 


0.8276 


0.8278 


0.826(2) 


0.18 


6S1/2 


5P1/2 


4.4346 


4.4565 


4.3672 


4.3504 


4.3606 


4.3529 


4.3529 


4.361(8) 


0.18 


6*1/2 


5p 3 /2 


6.3311 


6.3601 


6.2388 


6.2144 


6.2267 


6.2177 


6.2159 


6.23(1) 


0.17 


6S1/2 


6fl/2 


10.7169 


10.6885 


10.4530 


10.4853 


10.4977 


10.5108 


10.4952 


10.50(1) 


0.12 


6S1/2 


6f 3 /2 


15.1300 


15.0906 


14.7554 


14.8016 


14.8201 


14.8379 


14.8166 


14.82(2) 


0.12 


7S1/2 


4P1/2 


0.3140 


0.3276 


0.3193 


0.3183 


0.3182 


0.3186 


0.3187 


0.3183(6) 


0.20 


7S1/2 


4p 3 /2 


0.4456 


0.4644 


0.4530 


0.4515 


0.4514 


0.4520 


0.4520 


0.4515(9) 


0.20 


'Sl/2 


5^1/2 


1.1176 


1.1301 


1.1057 


1.1056 


1.1061 


1.1075 


1.1073 


1.106(1) 


0.12 


7S1/2 


5p 3 /2 


1.5845 


1.6020 


1.5676 


1.5675 


1.5681 


1.5702 


1.5699 


1.568(2) 


0.13 


7S1/2 


6P1/2 


7.4782 


7.4912 


7.3825 


7.3492 


7.3629 


7.3587 


7.3582 


7.363(5) 


0.06 


7S1/2 


6p 3 /2 


10.6693 


10.6877 


10.5391 


10.4913 


10.5071 


10.5040 


10.5009 


10.507(6) 


0.06 


7S1/2 


7P1/2 


16.0333 


16.0184 


15.7262 


15.7723 


15.7894 


15.8043 


15.7856 


15.79(1) 


0.09 


'I „ 

'Sl/2 


7p 3 /2 


22.6281 


22.6084 


22.1909 


22.2570 


22.2812 


22.3027 


22.2758 


22.28(2) 


0.10 


8S1/2 


5pi/ 2 


0.5862 


0.5941 


0.5827 


0.5829 


0.5830 


0.5840 


0.5839 


0.583(1) 


0.17 


8S1/2 


5p 3 /2 


0.8299 


0.8411 


0.8248 


0.8252 


0.8253 


0.8267 


0.8265 


0.825(1) 


0.17 


8S1/2 


6P1/2 


1.7969 


1.8017 


1.7732 


1.7738 


1.7743 


1.7770 


1.7767 


1.774(3) 


0.15 


8S1/2 


6p 3 /2 


2.5448 


2.5547 


2.5112 


2.5123 


2.5129 


2.5167 


2.5163 


2.513(4) 


0.15 


8S1/2 


7P1/2 


11.2499 


11.2611 


11.1339 


11.0842 


11.0946 


11.0951 


11.0879 


11.095(7) 


0.06 


8S1/2 


7p 3 /2 


16.0443 


16.0563 


15.8877 


15.8166 


15.8313 


15.8309 


15.8219 


15.831(9) 


0.06 


8S1/2 


8P1/2 


22.3941 


22.3854 


22.0348 


22.0969 




22.1359 




22.14(4) 


0.18 


8S1/2 


8p 3 / 2 


31.5981 


31.5878 


31.0854 


31.1745 




31.2302 




31.23(6) 


0.18 


3d 3/2 


4P1/2 


3.0825 


2.9296 


2.2998 


2.4173 


2.4636 


2.4677 


2.4503 


2.464(13) 


0.54 


3fl3/2 


4p 3 /2 


1.3764 


1.3088 


1.0260 


1.0788 


1.0996 


1.1014 


1.0937 


1.100(6) 


0.54 


3d 3/2 


4/5/2 


2.6059 


2.5228 


1.6763 


1.8660 


1.9265 


1.9265 


1.9051 


1.927(21) 


1.11 


3c<3/2 




1.5216 


1.4530 


1.1073 


1.1655 


1.1917 


1.1937 


1.1846 


1.192(7) 


0.59 


3d 5 / 2 


4p 3 / 2 


4.1348 


3.9311 


3.0882 


3.2452 


3.3063 


3.3127 


3.2884 


3.306(18) 


0.54 


O J 

3£t5/2 


4/5/2 


0.6976 


0.6751 


0.4502 


0.5005 


0.5163 


0.5166 


0.5106 


0.516(6) 


1.10 


3^5/2 




3.1201 


3.0192 


2.0134 


2.2382 


2.3090 


2.3102 


2.2835 


2.309(25) 


1.10 


3^5/2 


5/5/2 


0.4072 


0.3888 


0.2970 


0.3124 


0.3192 


0.3199 


0.3174 


0.319(2) 


0.59 




5/7/2 


1.8214 


1.7388 


1.3284 


1.3972 


1.4278 


1.4308 


1.4193 


1.428(8) 


0.59 


4d 3/2 


4P1/2 


4.2159 


4.1495 


4.4121 


4.2636 


4.2818 


4.2565 


4.2821 


4.28(3) 


0.59 


4fl3/2 


A ^, 

4p 3 /2 


1.8990 


1.8686 


1.9867 


1.9203 


1.9280 


1.9171 


1.9281 


1.93(1) 


0.57 


4d 3/2 


4/5/2 


11.9755 


11.8972 


11.3300 


11.3461 


11.3552 


11.4033 


11.3479 


11.36(5) 


0.42 


4d 3/2 


5pi/ 2 


8.0685 


8.0369 


7.2981 


7.4342 


7.4325 


7.4916 


7.4704 


7.43(6) 


0.79 


4d 3/2 


5p 3 /2 


3.6014 


3.5878 


3.2554 


3.3168 


3.3162 


3.3426 


3.3320 


3.32(3) 


0.79 


4d 5 /2 


4p 3 / 2 


5.6913 


5.6015 


5.9527 


5.7547 


5.7786 


5.7451 


5.7791 


5.78(3) 


0.58 


A J 

4d 5 /2 


5p 3 /2 


10.8150 


10.7724 


9.7806 


9.9632 


9.9600 


10.0403 


10.0091 


9.96(8) 


0.81 


4d 5/2 


4/5/2 


3.2020 


3.1809 


3.0300 


3.0342 


3.0364 


3.0494 


3.0345 


3.04(1) 


0.43 


4d 5 /2 


A £ 

4/7/2 


14.3198 


14.2254 


13.5508 


13.5696 


13.5793 


13.6376 


13.5707 


13.58(6) 


0.43 


5d 3 /2 


5P1/2 


7.1941 


7.1786 


7.6248 


7.4598 


7.4916 


7.4316 


7.4808 


7.49(6) 


0.80 


5d 3/2 


5p 3 /2 


3.2436 


3.2357 


3.4371 


3.3630 


3.3763 


3.3503 


3.3715 


3.38(3) 


0.77 


5d 3/2 


4/5/2 


7.9026 


7.9162 


8.7463 


8.5477 


8.5564 


8.4871 


8.5512 


8.56(7) 


0.81 


5d 5 / 2 


5p 3 /2 


9.7176 


9.6962 


10.2951 


10.0746 


10.1162 


10.0367 


10.1018 


10.12(8) 


0.79 


5d 5 / 2 


4/5/2 


2.1094 


2.1135 


2.3342 


2.2815 


2.2841 


2.2653 


2.2828 


2.28(2) 


0.82 


5d 5 / 2 


4/7/2 


9.4333 


9.4515 


10.4388 


10.2029 


10.2149 


10.1307 


10.2088 


10.21(8) 


0.82 


5d 3 / 2 


5/5/2 


20.4001 


20.3664 


19.0555 


19.1963 


19.2049 


19.3158 


19.1938 


19.2(1) 


0.58 


54/2 


5/5/2 


5.4558 


5.4462 


5.0977 


5.1350 


5.1367 


5.1668 


5.1337 


5.14(3) 


0.59 


54/2 


5/7/2 


24.3991 


24.3564 


22.7979 


22.9645 


22.9721 


23.1067 


22.9588 


23.0(1) 


0.59 


6d 3/2 


5P1/2 


2.3088 


2.2923 


2.2439 


2.2293 


2.2284 


2.2369 


2.2368 


2.228(9) 


0.39 


6d 3/2 


5p 3 /2 


1.0338 


1.0263 


1.0032 


0.9969 


0.9965 


1.0005 


1.0004 


0.996(4) 


0.40 


6d 5 /2 


5p 3 /2 


3.1018 


3.0793 


3.0112 


2.9922 


2.9908 


3.0027 


3.0025 


2.99(1) 


0.40 
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TABLE III: Comparison of the present values of El dipole 
matrix elements with RCC calculations of Ref. [22] . The un- 
certainties in our values represent our best estimate of all pos- 
sible sources of uncertainties, i.e. they give estimated bound- 
ary values of these recommended results. The uncertainties 
in Ref. [22I ] values are numerical uncertainties resulting from 
the use of incomplete basis sets. Absolute values in atomic 
units are given. 



Transition 


Present 


STOs 


Ref. [22J 
GTOs 


Final 


4f>l/2 


-4s 


2.898(13) 


2.86 


2.90 


2.88(1) 


4p 3 /2 


-4s 


4.099(18) 


4.02 


4.09 


4.03(1) 


#1/2 


— 3d 3 / 2 


2.464(13) 


2.50 


2.41 


2.40(2) 


4p 3 /2 


— 3d 3 / 2 


1.100(6) 


1.12 


1.09 


1.09(1) 


4f>3/2 


— 3d 5 / 2 


3.306(18) 


3.36 


3.28 


3.22(4) 



the valence triple excitations to the wave function given 
by Eq. ©. 

We find that only two terms give dominant contribu- 
tions for all matrix elements considered in this work: 



Z ' — ^ ] {ZamPwmva + z maP vmwa ) 



^ ^ ^ (z-wmprnv "T" Zmv Prnw) ' 



(3) 



(4) 



where p mna b = Pmnab - Pnmab and z wv are lowest-order 
matrix elements of the electric-dipole operator. For most 
of the transitions considered in this work, term is 
the dominant term. In many cases, it is overwhelmingly 
dominant (by a factor of 3 or more) . To evaluate missing 
corrections to this term, we need to improve the values of 
the valence single-excitation coefficients p mv (66|. These 
excitation coefficients are closely related to the correla- 
tion energy SE V . If we introduce the self-energy operator 
S mt , (also referred to as correlation potential in some 
works) as 



^mv \£v ^m) Pri 



(5) 



then the correlation energy would correspond to the di- 
agonal term T, vv [67j |. Therefore, the omitted correla- 
tion correction can be estimated by adjusting the single- 
excitation coefficients p mv to the experimentally known 
value of the valence correlation energy, and then re- 
calculating the matrix elements using Eq. ([T]) with the 
modified coefficients 166] 



'SEl heoly ' 



(6) 



The 5E° xpt is defined as the experimental energy [56| 
minus the lowest order DF energy e v . We note that it 
is a rather complicated procedure that involves complete 
recalculation of the matrix elements with new values of 
the valence excitation coefficients. The scaling factors 



depend on the correlation energy given by the particular 
calculation. Therefore, the scaling factors are different 
for the SD and SDpT calculations, and these values have 
to be scaled separately. Generally, scaled SD and SDpT 
values are close together, as expected. The corresponding 
results are listed in Table [TT| with subscript "sc" . 

The term is not corrected by the scaling pro- 

cedure. However, it is dominant for very few transi- 
tions that give significant contributions to the atomic 
properties considered in this work. Essentially, the only 
large matrix elements where term Z^ is not dominant 
are 4s — Ap and 5s — 5p. In both of these cases, term 
Z^> is still of the same order magnitude as the term 
Z( a \ Therefore, we can establish the recommended 
set of values and their uncertainties based on the ratio 
R = Z^/Z^. We take the final value to be SD scaled 
if R > 1. Otherwise, we use SD as the final value. If 
0.5 < R < 1.5, we evaluate the uncertainty in term Z^ 
as the maximum difference of the final value and the 
other three all-order values from the SD, SDpT, SDsc, 
and SDpTsc set. Then, we assume that the uncertainty 
of all the other terms does not exceed this value and 
add two uncertainties in quadrature. If 1.5 < R < 3, we 
evaluate the final uncertainty as the max(SDsc-SD, SDsc- 
SDpT, SDsc-SDpTsc). If the term Z^ strongly domi- 
nates and R > 3, we evaluate the final uncertainty as 
max(SDsc-SDpT, SDsc-SDpTsc). We note that we have 
conducted numerous comparisons of all available data on 
various properties of many different monovalent systems 
with different t ype s of exp eriments in many other works 
(see ^illliSlillsllli-lzl and references therein) 
and found that such procedures do not underestimate the 
uncertainties. If fact, they may somewhat overestimate 
the uncertainties in some cases. 

The last column of Table [TTI gives relative uncertainties 
of the final values Z in %. We find that the uncertain- 
ties are 0.2-0.5% for most of the transitions. Larger un- 
certainties (0.8%) occur for some of the transitions such 
as 5dj — Afji . Our final results and their uncertainties are 
used to calculate the recommended values of the transi- 
tion rates, oscillator strengths, lifetimes, and polarizabil- 
ities as well as evaluate the uncertainties of these results. 

Two most recent calculations of the El matrix ele- 
ments between the low-lying states were carried out by 
Sahoo et al. jUJ using the relativistic coupled-cluster 
method (RCC) and by Mitroy and Zhang _23] using 
non-relativistic configuration interaction with a semi- 



22] 



empirical core potential (CICP) approach. Ref. 
includes comparison with earlier MBPT calculations 
[H 51]. Since [H |H[ only include low-order MBPT 
corrections, these calculations are substantially less com- 
plete than all-order coupled-cluster method used in our 
work and Ref. [22(. Therefore, we focus our discussion 
on the comparison of the present results with those of 
[22| . Since Ref. [23| presents non-relativistic calculations 
and lists oscillator strengths rather than matrix elements, 
we compare their results with our j-averaged oscillator 
strengths in the next section. We note that we use the 
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TABLE IV: Comparison of the ratios R = S{ ld% of the squares of the El matrix elements calculated in the present work in 
several approximations with results of Ref. [22 ]. 



di 




di 


DF 


MBPT2 


Present work 
MBPT3 SD SDsc 


SDpT 


SDpTsc 


STOs 


Ref. [22] 
GTOs Final 


4s - 4p 3/2 
4pi/2 - 3d 3/2 

4P3/2 ~ 3d 5/2 


4s - 
4p 3 /2 

4P3/2 


- 3d 3/2 

- 3d 3 /2 


2.000 
5.02 
9.02 


2.001 
5.01 
9.02 


2.001 2.001 2.001 

5.02 5.02 5.02 
9.06 9.05 9.04 


2.001 
5.02 
9.05 


2.001 
5.02 
9.04 


1.976 
4.98 
9.00 


1.989 1.958(17) 
4.89 4.85(12) 
9.06 8.73(27) 



same method as [2l| and our results for the transitions 
listed in [2l[ are the same. Therefore, we do not include 
separate comparison with that work. 

Our final values are compared with RCC calculations 
of Sahoo et al. jH in Tabic HTD We discuss this compar- 
ison in significant detail since both calculations are car- 
ried out using the couple-cluster method but differ signif- 
icantly in its implementation. All El transitions listed in 
[22| are included. The results of [22| listed in columns la- 
beled "STOs" and "GTOs" are obtained by two different 
calculations, one with Slater-type orbitals and another 
with Gaussian-type orbitals, respectively. In both cases, 
the number of partial waves was restricted to l max = 4, 
i.e. only s,p, d, /, and g orbitals were included. Table [TTT1 
illustrates significant basis set dependence in the results 
of Sahoo et al. jUj, 1.4-1.7% for the 4s - Ap transitions 
and 2.4-3.8% for the 3d — Ap transitions. Moreover, dif- 
ferent fine-stricture components have different basis set 
dependencies. It is not clear how that may be possible 
unless some additional basis set optimization was car- 
ried out in different way for all transitions. The final 
results from [22| include corrections from higher symme- 
try orbitals carried out using MBPT(2); the changes be- 
tween the GTOs/STOs values and final recommended 
results range from 0.7% to 4.3%. Unfortunately, the 
Ref. [22j does not explicitly state what terms are ac- 
counted for by MBPT(2). In the standard formulation 
of the perturbation theory [frij , second-order perturba- 
tion theory contains only random-phase approximation 
terms which are relatively small for the 3d — Ap transi- 
tions. The main contribution of the higher partial waves 
only appears starting from the third order and comes 
from so-called Brueckner-orbital terms [64|. Therefore, 
using the MBPT(2) to evaluate higher symmetry contri- 
butions should severely underestimate these terms. The 
uncertainties of the final values from [22J are numerical 
uncertainties that are estimated from higher symmetry 
orbital corrections and consistency of results carried out 
with different basis sets. They range from 0.2% to 1.4%. 
We note that these are only numerical uncertainties and 
do not include estimates of missing theory (such as other 
triple and higher-excitation contributions). The uncer- 
tainties in our values represent our best estimate of all 
possible sources of uncertainties, i.e. they give estimated 
boundary values of the recommended results. We refer 
the reader to recent review [7o| for further discussion of 
the differences between numerical and complete uncer- 
tainties. 



The same very large basis is used in all calculations 
carried out in this work. We use 70 basis set functions for 
all partial waves with l max < 6. Use of such a large basis 
set results in negligible numerical errors in our values. 
The contribution of the / > 6 partial waves to the 4s — Ap 
transitions is expected to be at the 0.05% level. The 
contribution of the I > 6 partial waves to the 3d — Ap 
transitions is accounted for by the scaling procedure. The 
correction is small since the entire scaling of the ab initio 
SDpT values (that also accounts for the corrections due 
to higher-order and non-linear excitations) is 0.7%. 

The other major differences between the present work 
and Ref. [22[ include treatment of non-linear contri- 
butions, triple excitations, and higher-excitation terms. 
The Breit correction calculated in [22j is negligible for 
the El matrix elements. Ref. [22[ include non-liner terms 
at the SD level. While we have not explicitly included 
non-linear terms in this work, they were estimated by 
adjustment of the correlation potential described above 
along with contributions from higher excitations. We 
have demonstrated in Ref. [66] that correcting correla- 
tion potential staring from either linearized single-double 
coupled-cluster (LCCSD) or CCSD approximation leads 
to the same results within the expected accuracy of the 
calculations. Our ab initio inclusion of the valence triple 
excitation (SDpT) is more complete than that of [22[ 
since we included triple corrections to both SE V and p mv 
equations while only SE V was corrected in (22j. We also 
estimated other higher-excitation corrections as describe 
above while no such estimates were done in Ref. 

Our values are in agreement with results of 
the Api/2 — As and 4p 3 / 2 — 3d^/2 transitions but disagree 
well beyond the uncertainties for the 4p 3 / 2 — 4s, Ap 1 ^ 2 — 
3c?3/ 2 , and Ap 3 / 2 ~ 3d 5 / 2 transitions. We find this to 
be rather irregular since the correlation corrections are 
known to contribute nearly the same relative amount for 
the different fine-structure transitions for such light ions 
as Ca + . 

To clarify this issue, we calculated the ratios R — d\jd\ 
of the squares of all three relevant pairs of matrix ele- 
ments in the lowest-, second-, and third-order of MBPT, 
and all all-order approximations used in this work. The 
values of the matrix elements used in calculating the ra- 
tios are listed in Table HH The comparison of all val- 
ues is given in Table IIVI As we expected, all ratios of 
our values calculated in all approximations, including the 
lowest-order DF values, are nearly identical to the non- 
relativistic values (2, 5, and 9, respectively) which are 
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TABLE V: Wavelengths A (A), transition rates A (s _1 ) and oscillator strengths (/) for transitions in Ca II calculated using our 
recommended values of reduced electric-dipole matrix elements Z flnal and their uncertainties. The relative uncertainties in the 
values of transition rates and oscillator strengths are the same. They are listed in column "Unc." in %. Numbers in brackets 
represent powers of 10. 



Transition 


A 


4 

A 


/ 




Unc. 


Transition 


A 


A 


/ 




Unc. 


4si/ 2 


4pi/ 2 


3969.6 


1.360[8J 


3.213 




0.90 


4d 3/2 


4/5/2 


8914.5 


6.146[7J 


1.098 


0] 


0.84 


4si/ 2 


4^3/2 


3934.8 


1.397[8J 


6.485 


-!] 


0.90 


4d 5/2 


4/5/2 


8929.8 


4.372[6J 


5.227 


-2] 


0.86 


5S1/2 


5P1/2 


11953.0 


2.302[7J 


4.931 


-1] 


0.32 


4d 5/2 


4/7/2 


8929.8 


6.559|7J 


1.045 


0] 


0.86 


5si/ 2 


5P3/2 


H842.2 


2.362[7] 


9.932 


1 1 

-f] 


0.34 


5^3/2 


5/5/2 


18824.7 


1.867[7J 


1.488 




1.16 


6S1/2 


6^1/2 


26265.3 


6. 161 [6J 


6.372 




0.24 


5(4/2 


5/5/2 


18855.6 


1.329[6J 


7.084 


-2] 


1.18 


6S1/2 


6P3/2 


26013.7 


6.320[6J 


1.282 


0] 


0.24 


5d 5 / 2 


5/7/2 


18855.6 


1.994[7J 


1.417 


nl 

0] 


1.18 


A 

4f>l/2 


OSl/2 


3707.1 


8. 550 [7] 


1.761 


-1] 


1.02 


4f>l/2 


4d 3/2 


3159.8 


2.944[8] 


8.812 


-1] 


1.18 


4^3/2 


0^1/2 


3738.0 


1.705 [8] 


1.786 


-1] 


0.92 


4^3/2 


A W 

4<l3/2 


3182.2 


5.843[7] 


8.871 


-2] 


1.14 


4pi/2 


DS1/2 


2198.5 


3. 226 [7] 


2.338 


-2] 


0.36 


4^3/2 


4c(5/2 


3180.3 


3. 506 [8] 


7.973 


-1] 


1.16 


A ^ 

4f>3/2 


DS1/2 


2209.3 


6.416[7] 


2.347 


-2] 


0.36 


A „ 

4^1/2 


5rf3/2 


2103.9 


8.102[7] 


1.075 


-1] 


2.58 


4pi/2 


'Sl/2 


1843.1 


1.639 [7] 


8.349 


-3] 


0.40 


4P3/2 


0(13/2 


2113.8 


1.592[7] 


1.067 


-2] 


2.62 


A 

4f>3/2 


'Sl/2 


1850.7 


3.258[7] 


8.365 


-3] 


0.40 


4f>3/2 


5rf 5 /2 


2113.4 


9. 580 [7] 


9.623 


-2] 


2.62 


4pi/2 


Q n 

OSl/2 


1691.8 


9. 528 [6] 


4.088 


-3] 


0.42 


5^1/2 


oa s /2 


8204.0 


5.149[7] 


1.039 


0] 


1.60 


A ^ 

4f>3/2 


OSl/2 


1698.2 


1.893 [7] 


4.092 


-3] 


0.42 


5P3/2 


5^3/2 


8257.0 


1.026 [7] 


1.048 


-1] 


1.54 


OPl/2 


DSi/2 


9857.5 


2.011[7] 


2.930 


-1] 


0.36 


5P3/2 


5a 5 /2 


8251.1 


6.152[7] 


9.419 


-1] 


1.58 


5P3/2 


DSi/2 


9934.1 


4.007[7] 


2.964 


-1] 


0.34 


5P1/2 


003/2 


5002.9 


2. 009 [7] 


1.508 


-1] 


0.78 


5P1/2 


7s 1/2 


5286.7 


8. 388 [6] 


3.515 


-2] 


0.24 


5P3/2 


003/2 


5022.5 


3.970[6] 


1.501 


-2] 


0.80 


5P3/2 


7Si/ 2 


5308.7 


1.665 [7] 


3.517 


-2] 


0.26 


5P3/2 


6^5/2 


5021.4 


2.386[7] 


1.353 


-1] 


0.80 


5P1/2 


8S1/2 


4207.4 


4.623[6] 


1.227 


-2] 


0.34 












5P3/2 


8S1/2 


A OOI O 

4221. o 


y.l lo[0\ 


1.225 


-2J 


U.34 


3^3/2 


4P1/2 


8664.0 


9.4o2[6J 


r nifl 

o.oiy 


-2J 


1 no 
1.08 


5P1/2 


9«l/2 


3740.4 


2. 859 [6] 


5.997 


-3] 


0.38 


3d 3 / 2 


4P3/2 


8500.4 


9.972[5] 


1.080 


-2] 


1.08 


5P3/2 


9«l/2 


3751.4 


5. 670 [6] 


5.982 


-3] 


0.38 


3d 5 /2 


4P3/2 


8544.4 


8.876[6] 


6.477 


-2] 


1.08 


5P1/2 


lOSl/2 


3486.6 


1.899 [6] 


3.462 


-3] 


0.54 


4d 3/2 


5^1/2 


27072.6 


2.820[6] 


1.550 


-1] 


1.58 


5P3/2 


10S1/2 


3496.2 


3. 768 [6] 


3.453 


-3] 


0.54 


4d 3/2 


5P3/2 


26510.9 


2.990[5] 


3.150 


-2] 


1.58 


6P1/2 


7Si /2 


20147.6 


6.715[6] 


4.087 


-1] 


0.12 


4d 5 /2 


5P3/2 


26646.6 


2.656[6] 


1.885 


-1] 


1.60 


6P3/2 


7Si /2 


20298.3 


1.337 [7] 


4.130 


-1] 


0.12 












6P1/2 


8Si/ 2 


10187.5 


3.016[6] 


4.693 


-2] 


0.30 


4/5/2 


5^3/2 


21434.8 


3. 766 [6] 


1.729 


-1] 


1.62 


6P3/2 


8S1/2 


10225.8 


5. 983 [6] 


4.689 


-2] 


0.30 


4/5/2 


5d 5 / 2 


21394.9 


1.799 [5] 


1.235 


-2] 


1.64 












4/7/2 


5d 5 /2 


21394.9 


3. 598 [6] 


1.852 


-1] 


1.64 



simply the ratios of the corresponding angular factors. 
The effect of the entire correlation correction to the ra- 
tio is negligible. We see no feasible explanation of the 
anomalous ratios in [22| and significant changes between 
their GTO's/STO's and final value ratios. In the case of 
the 4s — Ap transitions, ratio of the final values from [22| 
is 3cr away from the NR value. The values of the other 
ratios are only slightly outside of the numerical error but 
it is not clear what could cause such changes in ratios 
from initial values to the final ones. The contributions 
from the higher symmetry orbitals can not change these 
ratios when the same basis set is used for the nl states 
with different j. 

In summary, we expect our results for all of the prop- 
erties listed in this work to be more accurate than that 



The transition rates are calculated using 



of Ref. [22[ based on the detailed analysis above. 



B. Transition rates and oscillator strengths 

We combine recommended NIST energies [5(| and our 
final values of the matrix elements listed in Table [TT| to 
calculate transition rates A and oscillator strengths /. 



A ab = 



2.02613 x 10 18 S 



A :i 



2j a + 1 



(7) 



where the wavelength A is in A and the line strength 
S = d 2 is in atomic units. 

Transition rates A (s _1 ) and oscillator strengths (/) 
for the 55 np — n's, np — n'd, and nd — n' f transitions in 
Ca II are summarized in Table |Vl Vacuum wavelengths 
obtained from NIST energies are also listed for reference. 
The relative uncertainties of the transition rates and os- 
cillator strengths are twice of the corresponding matrix 
element uncertainties since these properties are propor- 
tional to the squares of the matrix elements. The un- 
certainties in per cent are listed in the column labeled 
"Unc." . 

The values of the j-averaged oscillator strengths ob- 
tained using our final values of the matrix elements and 
NIST energies are compared with theoretical results from 
Refs. [H, in Table Ell The values of the Ref. [U 
are obtained with semi-empirical approach that uses ex- 
perimental energy levels and experimental or theoretical 
core polarizabilities as an input and approximates the 
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TABLE VI: Comparison of the j-avera ged oscillator strengths 
with theoretical results from Refs. [231 T4lT] 



Transitions 


Present 


Ref. [41J 


Ref. [231 


4s — 4p 


0.970(9) 


0.9523 


0.9606 


4s — 5p 


1.25F-3] 


3 31-41 


1.72 [-3] 


5s — 5p 


1.486(5) 


1.4736 




4p — 5s 


0.178(2) 


0.168 




4p - 4d 


0.884(10) 


0.8741 


0.8685 


4p - 5d 


0.107(3) 


0.1076 




3d-4p 


0.0645(7) 


0.0572 


0.0660 


3d — 5p 


5.1(1.8)1-4] 


6.9[-4] 


3.8[-4] 


4d — 5p 


0.188(3) 


0.1865 




3d -4/ 


0.154(3) 




0.1599 


4d-4/ 


1.098(9) 







core potential by the Hartree- Slater method. In recent 
work, Mitroy and Zhang [23] used non-relativistic config- 
uration interaction with a semi-empirical core potential 
(CICP) approach. The CICP values are in good agree- 
ment with our results taking into account the accuracy of 
both calculations. Earlier and significantly less sophisti- 
cated semi-empirical calculations of [4l| appear to be less 
accurate as expected. 



Lifetimes and branching ratios 



We calculated lifetimes of the 5s, 6s, 7s, 4pj, 5pj, 6pj, 
4dj, 5dj, and 4/j states in Ca + using out final values of 
the dipole matrix elements and NIST energies [56]. The 
uncertainties in the lifetime values are obtained from the 
uncertainties in the matrix elements listed in Table [TXJ 
The present values are compared with available experi- 
mental [H HI iH, HU and theoretical [H HJ results in 
Table Ell] 

The value of branching ratio of the 4p 3 / 2 decay 
of Ca + calculated from transition rates given in Ta- 
bleEM(4p 3/2 -4s)/XV A(4p 3/2 -3d J -)=14.15(20) agrees 
with 2008 measurement, 14.31(5), reported in Ref. [38[ 
within our uncertainty. As we noted above, our val- 
ues for the 4s — 4p and 4p — 3d 3 / 2 matrix elements 
are the same as in Ref. [21| since the same method 
is used. Therefore, the agreement of our values for 
the three branching fractions measured in (38j j re- 
mains the same as listed in the experimental work [38| : 
R(4p 3/2 - 4s) = 0.9347(3) expt vs. 0.9340 th , R(ip 3/2 - 
3d 3/2 ) = 0.00661 (4)°^ vs. 0.00667 th , i?(4p 3/2 -3d 3/2 ) = 
0.00587(2) ox P t vs. 0.00593 th . The uncertainties in our 
values of the transition rates are about 1% for all three 
transitions. Therefore, the agreement of our central val- 
ues with experiment is significantly better than expected 
from our uncertainty estimates (the uncertainty in the 
ratio is about twice that of the uncertainties in the indi- 
vidual transition rates). In fact, the 4s — 4p 3 / 2 branching 
fraction agrees with experiment to 0.07% making sub- 
stantial (3%!) disagreement of our 4p lifetime values with 



TABLE VII: Comparison of the lifetimes (in nsec) of nlj states 
with other theory and experiment. Uncertainties are given in 
parenthesis. References are given in square brackets. 



Level 



Present 



Expt. 



Theory 



5si/ 2 
6si/ 2 

7Si /2 

4pi/2 
4P3/2 
5P1/2 
5P3/2 
6P1/2 
6P3/2 

4d 



3/2 
5/2 



4(1 

5^3/2 
54/2 

4/5/2 
4/7/2 



3.91(4) 
6.39(2) 
10.63(3) 

6.88(6) 

6.69(6) 

35.4(7) 

34.8(7) 

89(2) 

90(2) 

2.83(3) 
2.85(3) 
6.16(13) 
6.21(14) 

3.55(7) 
3.54(7) 



4.3(4) [26] 



7.098(20) [31] 
6.924(19) [31] 



2.9(3) [27] 
3.1(2) [26] 
4-3(2) [26] 



4.153 [41] 
6.766 [41] 
11.262 [41] 

6.978(56) [22] 
6.926(36) [22] 
36.200 [41] 
35.249 [41] 
100.254 [41] 
99.675 [41] 

2.868 [41] 
2.886 [41] 
6.148 [41] 
6.199 [41] 

3.895 [41] 
3.897 [41] 



1993 experiment that lists 0.3% accuracy even more puz- 
zling. Our calculation of the 4p lifetimes in K [34j agrees 
with experimental values [77jj to 0.13%. Moreover, our 
primary npj lifetime values agree with recent experiments 
for all other alkali [H, Sr+ (36|, and Ba + |37|. Our values 
for the ground state polarizabilities that are completely 
dominated by the primary ns — np matrix element con- 
tributions, agree with all recent experiments in Li, Na, 
Cs [761, Mg+ [H|, Si 3+ , and Ba+(33|- This issue already 
have been discussed in detail in both 



21| and [3J. It 



would be very interesting to see new measurement of the 
Apj lifetimes, 4s — Apj transition rates, ground state po- 
larizability, or other properties that allow to infer 4s— Apj 
matrix elements in Ca + . 

All other experimental values listed in Table IVIII are 
much older (1970-1975) measurements with low preci- 
sion. The values of the metastable 3dj state lifetimes 
calculated with our approach agree within the uncertain- 
ties with the recent experimental values [llj]. This cal- 
culation and comparison with experiment was already 
discussed in detail in [llj], and we do not repeat it here. 
We note that 1% RCC theoretical value for the 3<i5/ 2 
lifetime, 1.110(9)s [2(J, calculated by the same group 
as work [22| that we discussed at length in the matrix 
element section disagrees with both our value and ex- 
periment by 6%. Nevertheless, their 3d 3 / 2 lifetime is in 
agrement with both our theoretical and experimental val- 
ues. This demonstrates another significant inconsistency 
of the approach used in [2^, HH in calculations of proper- 
ties of the levels from the same fine-structure multiplet. 
The ratio of these lifetimes is affected very weakly by the 
correlation as discussed in detail in[l8j]. Our value of this 
ratio is 1.0259(9) 0, while the ratio between Ref. [H 
3d lifetimes is 1.068. We note that lowest-order DF ratio 
of these lifetimes is 1.0245. Therefore, no difference in 
the treatment on the correlation correction can explain 
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such anomalous ratio of these lifetimes. 



IV. STATIC GROUND-STATE MULTIPOLE 
POLARIZABILITIES OF CA II 

The static multipole polarizability a Ek of Ca + in its 4s 
ground state can be separated to a valence polarizability 
and a polarizability of an ionic core. For the 4s state, 
the dominant valence contribution is calculated using the 
sum-over-state approach 



, Ek 



1 



2k 



nlj 



\{nl 3 \\r k C kq \\±s)\< 



En 



E A 



(8) 



where Ckq (f ) is a normalized spherical harmonic and nlj 
is npj, ndj, and nfj for k = 1, 2, and 3, respectively [35| . 
The E2 and E3 matrix elements and their uncertainties 
are calculated following the same approach that we used 
in calculating electric-dipole matrix elements (see Sec- 
tion inrs}. 

Contributions to the ground-state state dipole, 
quadrupole, and octupole polarizabilities are presented in 
Table IvTlTI Dominant contributions are listed separately. 
The remainders of the sums are listed together. For ex- 
ample, row labeled "nd 3 / 2 " gives the combined contribu- 
tion of all nd^/2 terms with n > 7. The first terms (4p, 
3d, and 4/, respectively) in the sum-over-states for a E1 , 
a E2 , and a E3 contribute 99.7%, 59%, and 79%, respec- 
tively, of the total valence polarizabilities. The rapid con- 
vergence of the sum over states for a E1 has been empha- 
sized in many publications (for example, Refs. pjil lilfy). 
The sums in Eq. ([5]) converge much slower for the E2 
and E3 polarizabilities. Therefore, accurate evaluation 
of a large number of terms in the sums ([5]) is needed for 
these states. We use NIST energies from [56| and our 
final recommended values of the matrix elements to eval- 
uate terms with n < 13. We use theoretical SD energies 
and matrix elements to evaluate terms with 13 < n < 26. 
The remaining contributions to a Ek from basis functions 
with 27 < n < 70 are evaluated in the DF approximation. 
These remainders are very small. Even in the case of the 
E3 polarizability, which is the slowest one to converge, 
the tail remainder with n > 26 contributes only 13 a.u. 
which is 0.14% of the total valence polarizability. 

The electric-dipole core polarizability is taken to be 
3.26(3) a.u. based on the comparison of the coupled- 
cluster and experimental values listed in the review [76| . 
This value is essentially the same as the random-phase 
approximation result of 3.25 a.u. Since this value is the 
polarizability of the ionic Ca 2+ core, we need to account 
for the presence of the valence electron by adding a term 
a vc which in this case is equal to half of the core polar- 
izability contribution from the excitation to the valence 
4s shell. In the cases of the E2 and E3 polarizabilities, 
we evaluate core contributions in the random-phase ap- 
proximation [79[ . The core polarizabilities are small in 
comparison with the valence ones and their uncertainties 



TABLE VIII: Contributions to dipole, a m , quadrupole, a E2 , 
and octupole, a ES , polarizabilities (a.u.) of the Ca + ground 
states. All values are in atomic units. The uncertainties are 
given in parenthesis. 



, El 



4z?i /o 
^■yi/ z 


24.4(2) 


3g?q /o 

%JXM o 1 2. 


203(2) 


4 ft; In 


3017(31) 


4^3/2 


48.4(4) 


4^3/2 


125(1) 


5/5/2 


469(4) 


npi/2 


0.08(2) 


5d 3 / 2 


10.6(2) 


6/5/2 


141(1) 




0.14(3) 


6d 3 /2 


2.7(1) 


7/5/2 


59.5(4) 


Core 




< a 3/2 




°/5/2 


■^n 7('\ \ 


OVc 


-0.12(1) 


nd 3/2 


4.0(2) 


nfs/2 


131(19) 


Total 


76.1(5) 


3ds/2 


304(3) 


4/7/2 


4023(41) 






4d 5 /2 


187(2) 


5/7/2 


625(5) 






5^5/2 


15.9(2) 


6/7/2 


188(1) 






6^5/2 


4.1(1) 


7/7/2 


79.3(5) 






7d 5 / 2 


1.6(0) 


8/7/2 


41.0(1) 






nd 5/2 


5.0(3) 


1/7/2 


175(24) 






Core 


6.9(3) 


Core 


34(4) 






Total 


871(4) 


Total 


9012(60) 


Ref. [23] 


75.49 


Ref. [23] 


875.1 


Ref. [23] 


8990 


Ref. [22] 


73.0(1.5) 











are negligible. We note that a vc terms are zero for the 
E2 and E3 polarizabilities since Ca + core contains no nd 
or nf states. 

Our final results for the ground-state multipole po- 
larizabilities are compared with other theoretical values 
[H HI in Table lyTITl The CICP values of Ref. [H 
are in remarkably good agreement with our results in all 
three cases. Interestingly, we differ by 3% in the main 
4s — 3dj contribution to the E2 polarizabilities (which 
is 523 a.u. in 23]), but agree in the final value. The 
difference between the present and Ref. [22[ El polariz- 
ability value results from the differences in the 4s — Apj 
matrix elements which we already discussed in detail in 
Section HIE1 



V. SCALAR AND TENSOR EXCITED STATE 
POLARIZABILITIES 

The valence scalar ao(t') and tensor 2 polarizabilities 
of Ca + in an excited state v are given by 



a (v) 



3(2j v + 1) 



E 

nlj 



\(v\\d\\nlj}\ 



E, 



nlj 



E„ 



(9) 



Q'2 



*£(-i) 

nlj 



i0j v (2j v - 1) 



3(> + l)(2^ + l)(2j;,+3) 
3 / jv 1 j \ \(v\\d\\nlj)\ 2 



1 in 2 J E n ij — E v 



(10) 



The ionic core polarizability discussed in the previous 
section has to be added to the valence term given by 
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TABLE IX: Contributions to the 5s, 6s, 7s, 8s, 4p,, 5pj, 3d,-, and Adj scalar polarizabilities of Ca II in ag. Uncertainties are 
given in parenthesis. The final results are compared with other theory [2^ . I23I ] . 



Contribution ao Contribution qq Contribution ao Contribution ao 





5s 






6s 




7s 




8s 




4pi/ 2 




-11.7(1) 


5pi/2 




-137.1(5) 


6P1/2 


-799(1) 


7Pl/2 


-3214(4) 


4p 3 /2 




-24.0(2) 


5P3/2 




-282(1) 
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-1639(2) 


7P3/2 


-6593(8) 
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339(1) 
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2118(5) 


7P1/2 


8894(21) 


8P1/2 


29153(125) 


5P3/2 




671(2) 


6P3/2 




4180(10) 


7P3/2 


17532(42) 


8P3/2 


57432(248) 


Other 




3.6 


Other 




3.8 


Other 


5 


Other 


-22 


Total 




978(3) 


Total 




5882(11) 


Total 


23990(50) 


Total 


76760(280) 


Ref. [23] 




983.5 
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4s 




-24.39(22 J 


OS 




-339(1] 


A ~> 

4pi/ 2 


-21.2(3] 


4P3/2 


or n / o\ 
-2o.y(3 ] 


5 s 




1 1 66(1 2] 


6s 




137 


u Kl /2 


547(9) 


5^3/2 


645(1 0) 


6s 




0.54 


7s 




5 


4P3/2 


-4.3 


«P3/2 





n s 




0.46(1) 


ns 




3 


5P3/2 


107(2) 


4/5/2 


20.1(2) 


3^3/2 




-38.47(41) 


4d 3/2 




-1094(17) 


np 





n/5/2 


0.1 


4d 3 /2 




42.38(50) 


5d 3/2 




337(5) 


4/5/2 


421(4) 


4/7/2 


402(3) 


5d3/2 




2.29(6) 


6c?3/2 




18 


0/5/2 


0.1 


5/7/2 


0.1 


nd 3/2 




1.51(8) 


nd 3 / 2 




11(1) 


n/5/2 


2.3(1) 


n/7/2 


2.2(1) 


Core 




3.26(3) 


Core 




3 


Core 


3.3 


Core 


3.3 


Total 




-0.75(70) 


Total 




-920(18) 


Total 


1054(10) 


Total 


1046(11) 


Ref. [23] 




-2.032 


Ref. [23] 




-1135 


Ref. [23] 


1209 


Ref. [23] 


1209 




4^3/2 






5P3/2 




3^3/2 




3^5/2 




A r, 
4S 




OA 10/ OOA 

-z4.1o(zz J 


OA 




-ooo(l ) 


4P1/2 


in OA /01 "S 
ly.z4(zl ) 


4P3/2 


OO 7Qi'0£.\ 

zz. /o(zo ) 


OS 




1 no f' 1 1 ^ 
lz.Uz(ll J 


OA 




141 


npi/2 


n no 
U.Uz 


np 3 /2 


n no 
U.Uo 


OS 




U.OO 


f s 




r 



4P3/2 


0. fb(4J 


A f 

4/5/2 


n 1 
U.lz 


ns 




0.46(1) 






3 


«P3/2 


0.01 


n/5/2 


0.17 


o«3/2 




-o. /0(4) 


4a 3 / 2 




1 A7/ 0\ 

-10Y(z ) 


A f 

4/5/2 


z.oU(O) 


4/7/2 


on/ c\ 

z.oy(o ) 


A 

4^3/2 




4.00(D) 


r J 

003/2 




34(1) 


0/5/2 


n 1 ( 1 \ 
U.ol(l ) 


5/7/2 


n 77/ 1 \ 

0.77(1) 


0«3/2 




U.ZO(± ) 


D«3/2 




9 
z 


DJ5/2 


n ^7? 1 ^ 

U.O / (1) 


0/7/2 


u.00^1 j 


nd 3/2 




0.15(1) 


nd z/2 




1 


7/5/2 


0.20(1) 


7/7/2 


0.19(1) 


3d 5 /2 




-34.17(37) 


4d 5/2 




-967(16) 


(8 - 12)/ 5 /2 


0.33(1) 


(8 - 12)/ 7/2 


0.31(1) 


4d 5 /2 




38.85(45) 


5^5/2 




309(5) 


(13 - 26)/ 5/2 


1.42(4) 


(13-26)/ r/2 


1.39(4) 


5^5/2 




2.07(5) 


6^5/2 




16 


nfs/2 


0.32(19) 


n/7/2 


0.27(15) 


nd 5/2 




1.22(6) 


nd 5/2 




9 


Core 


3.26(3) 


Core 


3.26(3) 


Core 




3.26(3) 


Core 




3 




-0.23(1) 




-0.23(1) 


Total 




1.02(64) 


Total 




-886(16) 


Total 


32.0(3) 


Total 


31.8(3) 


Ref. [23] 




-2.032 


Ref. [23] 




-1135 


Ref. [23] 


32.73 


Ref. [23] 


32.73 














Ref. [22] 


28.5(1.0) 


Ref. [22J 


29.5(1.0) 



Eq. ((9]) and corrected for the presence of the respective 
valence electron (term a vc ). This core correction a vc 
term is negligible for all excited states that we considered 
with the exception of the 3c? 3 / 2 and 3d 5 / 2 states. It is 
calculated in the RPA. 

The dipole polarizability calculations are carried out in 
the same way as the calculations of the multipole polar- 
izabilities discussed in the previous section. We list the 
contributions to the 5s, 6s, 7s, 8s, 4pj, 5pj, 3d,, and Adj 
scalar polarizabilities of Ca + in Table IIX1 The dominant 
contributions are listed separately. The remaining contri- 
butions are grouped together. For example, ll nd s / 2 " con- 
tribution includes all of the ndz/2 terms excluding only 
the terms that were already listed separately. 

The Tabic IIXI illustrates very fast convergence of the 
ns level polarizabilities which are dominated by the cor- 



responding np contributions. The (n — l)p term contribu- 
tions are significant, while all of the other contributions 
are very small owing to fast convergence of the sums. 
We use our recommended values for the 7p, and 8p en- 
ergy levels, for which we did not find the experimental 
values. The uncertainties in these recommended energy 
values are included when the polarizability uncertainties 
are calculated. The uncertainties of the final polarizabil- 
ity values are obtained by adding the uncertainties of the 
individual terms in quadarture. 

When calculating polarizabilities of the np and nd 
states, we calculated the terms with n < 26 using the 
all-order approach. The terms with n > 10 are calcu- 
lated using our calculated recommended values of the El 
matrix elements and the experimental energies [56j . The 
terms with 10 < n < 27 are calculated using SD energies 
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TABLE X: Contributions to the 4p 3 / 2 , 5p3/2, 3d-,, and Adj 
tensor polarizabilities of Ca II in ajj. Uncertainties are given 
in parenthesis. The final results are compared with other 
theory [Hill. 



Contribution 
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Ref [221 
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-22 45C5) 


Ref. [23] 


-17.64 


Ref. [23] 


-25.20 




4d 3 /2 




4rf 5 /2 


4pi/ 2 


21.2(3) 


4p 3 /2 


25.9(3) 


5pi/ 2 


-547(9) 


5P3/2 


-645(10) 


4p 3 /2 


-3.46(4) 


"P3/2 


-0.02 


5P3/2 


85(1) 


4/5/2 


23.0(2) 


np 


-0.01 


nfs/2 


0.13 


4/5/2 


-84.1(7) 


4/7/2 


-143(1) 


5/5/2 


-0.03(1) 


5/7/2 


-0.05(2) 




-0.46 


1/7/2 


-0.79 


Total 


-529(9) 


Total 


-740(10) 


Ref. [23] 


-615.9 


Ref. [23] 


-879.8 



and matrix elements. The remainders are evaluated in 
the DF or RPA approximations. The remainders in the 
3d — nf sums are treated more accurately as described 
below. 

We find that the scalar polarizabilities of the 4p 3 / 2 and 
4pi/2 states are anomalously small owing to a very pre- 
cise cancelations of the various contributions. This fact 
was already pointed out by Mitroy and Zhang 23j. Our 



uncertainties of these polarizability values are very large 
because of these severe cancelations. Such cancelations 
are not observed for the bp polarizabilities, where 5p — 4d 
contributions strongly dominate. 

The case of the 3d 5 / 2 polarizability is particularly in- 
teresting owing to its importance for the calculation of 



the blackbody radiation shift in the optical frequency 
standard with Ca + ion. Ref. [2l| points out that the 
sum over the n/7/2 states converges very slowly making 
accurate calculation of these contributions difficult. In 
this work, we have explored several different approaches 
to the accurate calculation of this sum and obtained con- 
sistent results in all cases. 

First, we calculated terms with n < 26 in the all- 
order approximation, and determined that these terms 
contribute 5.41(16) a.u. We find that even with so many 
terms included, the remainder is still 0.43 a.u. in DF ap- 
proximation and 0.32 a.u. in the RPA, which is a signif- 
icant fraction of the total n/7/2 contribution. Therefore, 
we estimate the accuracy of the DF approximation by cal- 
culating the main terms with n < 26 in the DF approxi- 
mation as well. We find that DF approximation overesti- 
mates the polarizability contributions from highly-exited 
state by about 38% and adjust the DF value accordingly. 
The entire adjustment is taken to be the uncertainly of 
the n > 26 remainder. Therefore, our final value for 
71 > 26 77/7/2 terms is 0.27(15), and the total 77/7/2 sum 
is equal to 5.67(22) a.u. 

Second, we carry out the calculation of the highly- 
excited states by a different approach to verify that the 
all-order calculation of such highly-excited states does 
not introduce unexpected errors. We compare the con- 
tributions with 6 < 77, < 13 calculated in DF, RPA, and 
all-order approximations. The wave functions of all these 
states fit inside of our 220 a.u. cavity, and the all-order 
method is definitely reliable for these states. We establish 
that DF overestimates the results by 40-44% and RPA 
overestimates the results by 37-41%. Then, we use these 
percentages to adjust the contributions for the 77 > 12 
states calculated in the DF and RPA approximations. 
Adding these adjusted remainders to the all-order terms 
with 77 < 13 gives 5.63 a.u. (DF) and 5.53 a.u. (RPA) 
for the total 77/7/2 sum. These values are consistent with 
our result 5.67(22) obtained above. Such accurate eval- 
uation of this sum allows us to reduce the uncertainty 
of the 3d 5 /2 polarizability by a factor of 3 in compari- 
son with the previous calculation of this quantity using 
the all-order sum-over-states approach [2l| . The calcula- 
tion of the 71/5/2 contribution to the 3(^/2 polarizability 
is carried out using the same method. In all the other 
cases, the contributions of the terms with 71 > 26 are 
very small in comparison with the other terms and DF 
approximation is sufficiently accurate. 

We compare our values of the scalar polarizabilities 
for the 3d 3/2 and 3d 5/2 with RCC || and CICP [H 
theoretical calculations in Table IIXI Ref. [22| estimates 
their numerical (basis set truncation) uncertainty to be 
3.5%. This estimate does not account for the uncertainty 
owing to the missing correlation correction. Our 3^ po- 



larizabilities differ from values of [22[ by 7% and 11%. It 
is well known (see [76| and references therein) that the 
polarizabilities are very sensitive to the problems with 
the basis set completeness in the coupled-cluster calcula- 
tions such as RCC calculation of Ref. [22[. We note that 
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Ref. HI] STO values (31.6 and 32.5) are very close to our 
results. 

Since the CICP calculation 23] is non-relativistic, we 
list their values for both fine-structure states. Our 5s 
and 3d results are in good agreement with CICP calcu- 
lation. The agreement is rather poor for the 4p, 4d, and 
5p states. It is expected for the 4p states owing to severe 
cancelations discussed above but somewhat surprising for 
the other two states. 

We list the contributions to the tensor polarizabilitics 
of the Ca + in 4p 3 / 2 , 5j>3/2> 3dj, and 4dj states in Ta- 
ble [Xj Tensor polarizability calculations are carried our 
in the same way as the scalar polarizability ones. The 
same designations are used in Table [X] as in the scalar 
polarizability Table IIX1 The final values are compared 
with RCC 22] and CICP 23J theoretical calculations. 
We multiply the 3d and Ad non-relativistic values of 23 1 
by 7/10 to compare these values to our 3ti 3 / 2 and 4d 3 / 2 
tensor polarizabilities (see [761 ] for explanation of this con- 
version factor). The differences between the present and 
[22l [23| theoretical tensor polarizability values are similar 
to those for the scalar polarizabilities for these states. 



VI. BLACKBODY RADIATION SHIFT IN CA+ 
OPTICAL FREQUENCY STANDARD 

The electrical field E radiated by a blackbody at tem- 
perature T, as given by Planck's law, 



TABLE XI: Hyperfine constants A (in MHz) in 43 Ca+ (I = 



E 2 {u)duj 



8a J 



3 dw 



7r exp(w//c.BT) — 1 ' 



(11) 



induces a nonresonant perturbation of the optical tran- 
sition at room temperature [80| . The frequency shift 
of an atomic state due to such an electrical field is re- 
lated to the static electric-dipole polarizability ao by (see 
Ref. M) 



Av = "(831.9 V/m) 2 ( T ^ K ^ 



\ 300 ) 



a (l + r)), (12) 



where r\ is a small dynamic correction due to the fre- 
quency distribution. Only the electric-dipole transition 
part of the contribution is considered in the formula 
above because the contributions from Ml and E2 transi- 
tions are suppressed by a factor of a 2 [Hj]. The overall 
BBR shift of the Ca + 4s — 3d 5 / 2 clock transition fre- 
quency is then calculated as the difference between the 
BBR shifts of the individual levels involved in the tran- 
sition: 

A B br(4s - 3d 5 / 2 ) = -^[a (3d 5/2 ) - a (4s 1/2 )] 

x (831.9V/m) 2 (^T) 4 .(13) 

The tensor part of polarizability is averaged out due 
to the isotropic nature of the electric field radiated 



7/2, » 


= -1.31727). The SD and SDpT all-order results 


are compared with theoretical 551] 


and experimental [5ll, |53|] 


values. 










Level 


A DF 




^SDpT 


Th [551 F.-srnt [51 5S1 


4s 


-587.39 


-818.82 


-801.28 


-805.348 -806.402072 


5s 


-183.99 


-239.19 


-236.30 




6s 


-81.25 


-103.29 


-102.33 




7s 


-42.91 


-53.96 


-53.53 




4 Pl/2 


-101.48 


-148.26 


-144.96 


-143.068 -145.4(1) 


5pi/2 


-37.43 


-51.28 


-50.71 




6P1/2 


-17.91 


-24.14 


-23.81 




TPl/2 


-9.94 


-13.25 


-13.08 




4P3/2 


-19.64 


-31.04 


-30.34 




5P3/2 


-7.25 


-10.69 


-10.58 




6P3/2 


-3.47 


-5.04 


-4.97 




TP3/2 


-1.93 


-2.76 


-2.73 




„ , 


t>t).iU 


4:0.00 


47 


47 894 47 ^C?} 


4^3/2 


-8.01 


-9.32 


-9.35 




5d 3 / 2 


-3.36 


-3.86 


-3.88 




6^3/2 


-1.75 


-1.99 






7d 3 /2 


-1.03 


-1.17 






3^5/2 


-14.14 


-4.71 


-4.24 




4d 5 / 2 


-3.41 


-3.12 


-3.06 




5d 5 / 2 


-1.43 


-1.51 


-1.49 




6d 5 / 2 


-0.748 


-0.831 






7d 5/2 


-0.441 


-0.502 






4/5/2 


-0.151 


-0.165 


-0.163 




5/5/2 


-0.079 


-0.089 






6/5/2 


-0.046 


-0.053 






7/5/2 


-0.029 


-0.034 






4/7/2 


-0.084 


-0.044 


-0.044 




5/7/2 


-0.044 


-0.011 






6/7/2 


-0.026 


-0.002 






7 /7/2 


-0.016 


0.000 







by the blackbody. 



4s and 3d 5 / 2 static 



Substituting out values for the 
polarizabilities into Eq. (|13[) . we 
obtain 0.3815(44) Hz for the BBR shift. We note 
that atomic units for a are converted to SI units via 
a/h[Rz/(V/m) 2 } = 2.48832 x 10- 8 a[a.u.], where the con- 
version coefficient is Aire^a^/h and Planck constant h is 
factored out. 

We estimate the dynamic corrections to be 77=0.0012 
and 71=0.0044 for the 4s and 3d 3 / 2 states, respectively, 
following Ref. [8l[. The resulting dynamic correction to 
the BBR shift is —0.0004 Hz and our final value is 

A BB r(4s - 3d 5/2 ) = 0.3811(44) Hz. 

The value is the same for different Ca + isotopes within 
its accuracy. The third-order independent polarizability 
of the ground state is evaluated in the last section of this 
paper. Its contribution is several orders of magnitudes 
smaller than the second-order value and can be omitted 
in evaluating BBR shift in the optical standard. 

The present value is consistent with other calculations, 
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TABLE XII: Hyperfine constants B (in MHz) in 43 Ca + . Nuclear quadrupole moment Q is taken to be equal to -0.044(9) barns 
(1 b=10" 24 cm 2 ) The SD and SDpT values are compared with theory [55|| and experiment [53l ]. 



Level 


S DF 

Q 


gsn 


gSDpT 

Q 






^SDpT 


^ [551 


B cxpt [53J 


4^3/2 


96.69 


153.99 


150.81 


-4.25 


-6.78 


-6.64 


-151.798 


-6.9(1.7) 


5P3/2 


35.58 


52.40 


51.82 


-1.57 


-2.31 


-2.28 






6^3/2 


16.93 


24.37 


24.06 


-0.75 


-1.07 


-1.06 






3d 3 /2 


54.24 


67.01 


65.24 


-2.39 


-2.95 


-2.87 


-68.067 


-3.7(1.9) 


4d 3 / 2 


12.93 


17.37 


17.24 


-0.57 


-0.76 


-0.76 






54/2 


5.31 


7.28 


7.24 


-0.23 


-0.32 


-0.32 






34/2 


76.86 


95.19 


92.69 


-3.38 


-4.19 


-4.08 


-100.208 


-3.9(6.0) 


4d 5/2 


18.33 


24.70 


24.51 


-0.81 


-1.09 


-1.08 






5d 5 /2 


7.54 


10.35 


10.30 


-0.33 


-0.46 


-0.45 







0.380(14) Hz [U, 0.37(1) Hz 22], and 0.368 Hz [Hj], but 
is three times more accurate. 



VII. HYPERFINE CONSTANTS FOR M CA+ 

Calculations of hyperfine constants are carried out us- 
ing the SD and SDpT all-order methods described in Sec- 
tion IIII Al A number of terms other than terms Z^ a > and 
Z(°) give significant contributions to the hyperfine con- 
stants. Therefore, scaling procedure described in Section 
IIII Al is not expected to produce more accurate values 
and is not carried out for the hyperfine constants. In 
Table IXI| we list hyperfine constants A for 43 Ca + and 
compare our values with available theoretical (EBJ and 
experimental data [H, [53[ . 

In this table, we present the lowest-order A° F , all- 
order ^4 SD , and ^4 sd p t values for the ns, np, nd, and nf 
levels up to n = 7. The magnetic moment of 43 Ca + used 
here (I = 7/2, fj, = -1.31727) is taken from 82]. Our 
SDpT results are in very good agreement with experi- 
mental results for the ns and npx/2 states when experi- 
mental uncertainties are taken into account. The contri- 
butions from valence triple excitations are large for the 
hyperfine constants and have to be included for an accu- 
rate calculation. 

Hyperfine constants B (in MHz) in 43 Ca + are given 
in Table IXII1 Nuclear quadrupole moment Q is taken to 
be equal to -0.044(9) barns (1 b=10" 24 cm 2 ) [IpL The 
SD and SDpT values are compared with theory [55| and 
experiment [53| . 



VIII. HYPERFINE- INDUCED TRANSITION 
POLARIZ ABILITY OF THE 43 CA+ GROUND 
STATE 

We now turn to the calculation of the quadratic Stark 
shift of the ground-state hyperfine interval (F — 4 — 
F = 3) in 43 Ca + . The quadratic Stark shift is closely 
related to the blackbody radiation shift in the microwave 
frequency standards discussed, for example, in Refs. [83l — 
HH . Our calculation follows the methodology outlined in 



those works. 

The dominant second-order contribution to the polar- 
izability cancels for the transition between the two hy- 
perfine components of the 4s state. Therefore, the Stark 
shift of the hyperfine interval is governed by the the third- 
order F-dependent polarizability a F (0) . The expression 
for the c4 3) (0) is ||: 

4 3) (0) = i v / (2/)(27+l)(2J + 2)| ^ f }x 

9ifin(-lf +I+h (2T + C + R), (14) 

where gi is the nuclear gyromagnetic ratio, fi n is the 
nuclear magneton equal to 0.3924658 in 43 Ca + , I = 7/2 
is the nuclear spin, and j v = 1/2 is the total angular 
momentum of the atomic ground state. The formulas 
for the F- independent terms T, C, and R are given in 
Ref. These terms are similar to the polarizability 

sum-over-state expression but are more complicated. 

First, we calculate these values in the DF approxima- 
tion (in atomic units): 

2T df = 2.0018 x 10~ 4 , C* DF = 3.9507 x 10~ 7 , 
i? DF = 3.8838 x 10~ 4 . (15) 

Since the value of C DF is smaller than the T DF and 
i? DF by three orders of magnitude, we do not recalculate 
the C term using the all-order method. 

The expression for R is similar to that for a E1 but 
contains diagonal hyperfine matrix element: 

(4s||r||4s) SDpT = 3.9629 x 10~ 7 a.u. 

We use our all-order recommended values for the re- 
duced electric-dipole matrix elements described in Sec- 
tion IIII Al and their uncertainties to calculate the main 
terms in the T and R sums. We refer to these values 
as the "best set" values. Available recommended NIST 
energies [56[ are used for nl = 4s— 10s, 4p—6p, and SD en- 
ergies are used for the other states up to n = 26. The sum 
of R terms with n < 26 is equal to R = 3.772(34) x 10~ 4 . 
The remainder of the R sum is evaluated in the DF ap- 
proximation, i?„>26 = 3-0 x 10~ 8 , and is less than 0.01%. 



16 



TABLE XIII: Contributions to the mp sums of term 2T, 
m = 4 — 26. The main contribution y^ 26 calculated in the 
DF approximation is given in the column labeled "Main DF " 
in 10 4 a.u. The final values of the main contributions to 
the mp sums are given in the column labeled "Main flnal " in 
10 4 a.u. Accumulated values are given for both DF and final 
results. The ratio of the final and DF values for the main 
terms is given in the fourth column in %. The relative tail 
contribution X^I=27 calculated in the DF approximation is 
given in the last column. 



mp 


Main DF 


Main flnal 


Dif.(%) 


Tail (%) 


An 


1.965 


2.060(12) 


4.6 


1.8 


- 1 


1.965 


2.079(12) 


5.5 


1.8 


(yd 


1.969 


2.089(12) 


5.7 


1.6 


' if 


1.971 


2.094(12) 


5.9 


1.5 


Ho 

LJJJ 


1.972 


2.097(12) 


5.9 


1.5 


Qr)f) 


1.973 


2.098(11) 


6.0 


1.4 


lOp 


1.974 


2.100(12) 


6.0 


1.4 


lip 


1.974 


2.101(12) 


6.0 


1.4 


I2p 


1.974 


2.102(13) 


6.1 


1.4 


13p 


1.975 


2.102(13) 


6.1 


1.4 


Up 


1.975 


2.102(13) 


6.1 


1.4 


15p 


1.975 


2.102(13) 


6.1 


1.4 


16p 


1.975 


2.102(13) 


6.1 


1.4 


17p 


1.975 


2.103(13) 


6.1 


1.3 


18p 


1.976 


2.104(13) 


6.1 


1.3 


19p 


1.976 


2.104(13) 


6.1 


1.3 


20p 


1.976 


2.104(13) 


6.1 


1.3 


21p 


1.979 


2.109(13) 


6.2 


1.2 


22p 


1.979 


2.109(13) 


6.2 


1.2 


23p 


1.983 


2.118(13) 


6.4 


0.9 


24p 


1.986 


2.124(13) 


6.5 


0.8 


25p 


1.988 


2.129(13) 


6.6 


0.7 


26p 


1.990 


2.129(13) 


6.5 


0.6 



Term T contains two sums, over ns and over mpj. We 
evaluate main contributions, that include n < 26 and 
m < 26 using all-order matrix elements and NIST or 
all-order energies as described above. We find that the 
remaining contributions with n > 26 and m > 26, are 
very small. Table IXIIII illustrates the size of the remain- 
ders and accuracy of the DF approximation. We break 
down each mp term as 

70s \ 




nip 



27s 



and list contributions to the mp sums of term T, m = 
4 — 26 in Table IXIIII Each mp term is given by 



2(i 



n— 5 



(<is\\D\\mpj)(mpj\\D\\ns)(ns\\T\\^s) 
(E m p — E^s) (E ns — Ej s ) 



(16) 



where At is an angular factor. The main contribution 
zCn=5[---] calculated in the DF approximation is given in 
the column labeled "Main DF " in 10 4 a.u. The final values 
of the main contributions to the mp sums are given in 
the column labeled "Main final " in 10 4 a.u. Accumulated 



values are given for both DF and final results to illustrate 
the convergence of the mp sum. The ratio of the final 
and DF values for the main terms is given in the fourth 

- ^70 

column in %. The relative tail contribution X)n=27[-"] 
calculated in the DF approximation is given in the last 
column. The remainder is 0.6% of the main term and 
is equal to 0.013 x 10~ 4 a.u. Our final value for this 
term is 2T final = 2.142(13) x 10~ 4 . Combining these 
contributions, we obtain 



2T 



final 



c 



DF 



R bna[ = 5.918(36) x 10~ 4 a.u. (17) 



The F-dependent factor in Eq. flU) is equal to 0.4609 
for F — 3 and -0.3585 for F = 4. Using these values and 
the result from Eq. (|17p . we obtain 



a hf (4s) = c4 3 l 4 (0) - o^i 3 (0) = -4.850(29) x 10" 
The Stark shift coefficient k defined as Av 



„(3) 



a.u. 



kE 2 



a F=4(0) — a Fi=3(0) • Converting from atomic 



„(3) 



units, we obtain 

k = -2.425(15)xl0" 4 ; 



6.03(4) xl0^Hz/(V/m) 



We note that the lowest-order DF value is fc( DF ' = 
6.00xl0~ 12 Hz/(V/m) 2 . While values of both R and T 
terms change with the inclusion of the correlation correc- 
tion, it essentially cancels when these terms are added. 
The relative blackbody radiative shift (3 is defined as 



P 



2 1 

15 um 



(air) r 4 a h f(4si/ 2 ) 



(18) 



where i/y is the 43 Ca + hyperfine (F = 3 — F = 4) split- 
ting equal to 3225.6082864(3) MHz [U and T is tem- 
perature taken to be 300 K. Using those factors and our 
value of % (4s) , we obtain 



P = -2.6696 x 10- L2 a h{ (4s) = 1.29(1) x 10 
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IX. CONCLUSION 

A systematic study of Ca + atomic properties is carried 
out using high-precision relativistic all-order method 
where all single, double, and partial triple excitations of 
the Dirac-Fock wave function are included to all orders 
of perturbation theory. Energies, El, E2, E3, matrix 
elements, transition rates, lifetimes, A and B hyperfine 
constants, El, E2, and E3 ground state polarizabilities, 
scalar El polarizabilities of the 5s, 6s, 7s, 8s, Apj, 5pj, 
3dj, Adj states, and tensor polarizabilities of the 4p 3 / 2 , 
5p 3 /2, 3dj, and 4dj states are calculated. We evaluate 
the uncertainties of our calculations for most of the 
values listed in this work. The blackbody radiation 
(BBR) shift of the 4s — 3(^/2 clock transition in Ca + 
is calculated to be 0.381(4) Hz at room temperature, 
T = 300 K improving its accuracy by a factor of 3. The 
quadratic Stark effect on hyperfine structure levels of 
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Ca + ground state is investigated. These calculations 
provide recommended values critically evaluated for 
their accuracy for a number of Ca + atomic properties 
useful for a variety of applications. 
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